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Abstract 

We consider perfect fluid bodies (“stars”) in general relativity that are axisymmetric, asymptot¬ 
ically flat, and that admit a maximal hypersurface. We show that configurations that extremize 
the total entropy at fixed ADM mass, ADM angular momentum, and total particle number are 
stationary with circular flow. For such stars, this establishes that thermodynamic equilibrium 
implies dynamic equilibrium. 
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I. INTRODUCTION 


In the description of self-gravitating flnid bodies—“stars”—in general relativity, the fnn- 
damental conserved qnantities that characterize the global state of a star can be taken to 
be the ADM mass M, the ADM angular momentum J, and the total number of particles 
N. The Einstein-perfect fluid equations predict an additional conserved global quantity: 
the total entropy, S. However, S is not on equal footing with M, J, and N. The per¬ 
fect fluid description is an idealized one in which dissipative processes—i.e., viscosity, heat 
conduction, and diffusion—cannot occur. In a more realistic fluid description that takes 
dissipative processes into account, a star that is not in equilibrium would be expected to 
have S increasing with time, while M, J and N remain constant. One would expect the 
star to eventually settle down to a stationary end-state that is a (local) maximum of S for 
the given values of M, J, and N. 

Suppose, now, that we are given initial conditions for a non-stationary perfect fluid star. 
Imagine that we “turn on” a small amount of dissipation at the initial time, wait a moment, 
and then “turn off” the dissipation. We would expect S to increase during the time that 
the dissipation is “on”. Thus, we will have produced a perfect fluid state “nearby” the 
initial state that has higher entropy. This suggests that extrema of S at fixed M, J, and N 
must be stationary. Furthermore, Lindblom [l| has shown that stationary dissipative fluid 
stars are necessarily axisymmetric with “circular flow,” meaning that the fluid 4-velocity 
is a linear combination of the stationary and axial Killing vector fields, and Thus, 
if we consider a stationary perfect fluid star that lacks axisymmetry or circular flow and 
we do our experiment of turning on some dissipation, we would again expect the entropy 
to increase. This line of reasoning suggests we make the following conjecture: Any perfect 
fluid star that is an extremum of entropy—i.e., hS* = 0 for all linearized solutions that have 
6M = 6 J = 6N = 0—is stationary and axisymmetric with circular flowll] 

In this paper we give a proof that any axisymmetric perfect fluid star that admits a max¬ 
imal (i.e., zero trace of the extrinsic curvature) hypersurface and that is at an extremum of 
entropy (at hxed M, J, and N) is stationary with circular flow. In other words, we will prove 
the above conjecture for stars that are assumed from the beginning to be axisymmetric, and 

^ One may furthermore conjecture that any extremum of entropy that is in fact a local maximum—i.e., 
6^S < 0 for all perturbations that keep M, J, and N fixed to both first and second order—is dynamically 
stable. In fact, under the assumptions of stationarity, axisymmetry, and circular flow, this was shown to 
be true in Q (and it was earlier shown in the case of static, spherically symmetric radiation stars by [^, 
and in the case of static, spherically symmetric barotropic stars by 
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that admit a maximal hypersurface. As we will discuss in section IIIIl the latter assumption 
is not believed to exclude any physically reasonable spacetimes. 

The question of proving the above conjecture has been previously addressed in the lit¬ 
erature in the case where spherical symmetry is assumed, in which case J is removed from 
the problem. In particular, Sorkin, Wald, and Zhang considered radiation fluid enclosed 
in a spherical box. They argued—but as we will describe below, did not prove—that any 
conhguration that admits a maximal hypersurface and that extremizes S at hxed M (there 
is no notion of N for a radiation fluid) is static^. Our stationarity proof is modeled after 
their staticity argument, so we will review it here. Afterwards, we will explain what the 
obstacles are in extending their method to the case of axisymmetric rotating stars, and how 
we overcome these obstacles. 

The basic idea behind the argument of Sorkin, Wald, and Zhang is as follows: If there is 
any bulk fluid flow, it should be possible to construct a perturbation that slows down this 
flow and converts its energy into thermal energy, thereby raising the entropy. Constructing 
perturbed initial data that does this is nontrivial because the constraint equations must 
be satished. They start by assuming the existence of a maximal hypersurface S. The 
momentum constraint on such a surface is 

= -SvrPi, (1) 


and the Hamiltonian constraint is 


R{h) = IQ'ke + KijK^^, (2) 

where hij is the spatial metric, Kij is the extrinsic curvature. Pi is the fluid momentum 

density, and £ is the fluid energy density (as viewed by an observer moving normal to S). 

They construct their perturbation on E by hrst doing a “momentum scaling” which scales 

down Kij and Pi by a constant factor while keeping hij hxed. As hij is unchanged, so is 

M. The momentum scaling doesn’t change the momentum constraint ([T]), but it decreases 

^ Other related results have been proven in the spherically symmetric case. In Q it was shown for con¬ 
figurations of a radiation fluid that admit a time-symmetric (i.e., zero extrinsic curvatu^ hypersurface 
that being an extremum of entropy is sufficient and necessary for being static, and Gao extended the 
“sufficient” part of this result to allow for a general 2-parameter equation of state (the “necessary” part 
holds only for 1-parameter equations of state). Fang and Gao Q obtained some related results without 
assuming spherical symmetry, but they assumed from the start that the spacetime metric is static and 
that the redshifted temperature is uniform. 
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the term KijK^^ in the Hamiltonian constraint ([2]). They compensate for this by increasing 
e, which can be shown to increase the entropy density of the radiation fluid. Thus, if the 
conhguration is an extremum of S at hxed M, we must have Kij = 0—i.e., E must be 
a time-symmetric hypersurface. Additionally, by ([I]), the fluid momentum density must 
be zero, implying that the fluid 4-velocity is normal to E. But existence of one maximal 
hypersurface E implies j?] a unique foliation by maximal hypersurfaces in a neighborhood of 
E, so it follows that a neighborhood of E can be foliated by time-symmetric hypersurfaces. 
But this means that the neighborhood, and therefore the entire spacetime, is static. 

Where their argument fails to be a proof is in the step where they compensate for the 
decrease in KijK^^ with an increase in £. Although this works to keep the Hamiltonian 
constraint satished inside the fluid body, it does not work outside the body where e = o|^ 

The basic idea of our stationarity proof in the axisymmetric case, in parallel with the 
basic idea of the above staticity argument, is as follows: If the fluid flow has any “non¬ 
circular” part—i.e., a part perpendicular to the axial Killing vector held 0“—it should be 
possible to construct a perturbation that slows down this non-circular how and converts its 
energy into thermal energy, thereby raising the entropy. Only the circular part of the how 
contributes to the star’s angular momentum, so such a perturbation should not change J. 

We will construct a perturbation analogous to the one of described above to use in the 
axisymmetric case with general equation of state to prove our stationarity theorem. There 
are three main obstacles to doing this: (i) In addition to M, the angular momentum, J, 
must now be held hxed under the perturbation, (ii) The particle number, N^ must also be 
held hxed. (hi) The perturbation must satisfy the constraint equations everywhere on E, 
including outside the star where £ = 0. 

Obstacle (ii) (keeping N hxed) is straightforward to overcome, since, as we will show, 
one has the freedom to keep the particle number density hxed while scaling down P, and 
increasing the entropy density. Overcoming obstacle (i) (ke( ' ^ ' ontrivial. 



however. We will do this—following methods developed 


vacuum 


spacetimes—by using the axisymmetry to decompose hij, Kij and Pi into “axial” and “polar” 


parts with respect to the axial Killing vector held. Our analog of the “momentum scaling” 
will be dehned to scale down the axial part of hij and the polar parts of Kij and P, by a 


^ Their argument can be made into a proof by “fixing up” their perturbation via a linearized conformal 
transformation of the initial data, as we will do below in the general axisymmetric case. 
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constant factor, while keeping the polar part of hij and the axial parts of Kij and P, fixed. 
As we will show, this keeps M and J nnchanged. 

To overcome obstacle (hi) we need a way to “fix np” changes in the Hamiltonian constraint 
other than by changing e to compensate. Following and 8], we will do so via a linearized 
conformal transformation of the initial data, i.e., a linearized version of the Lichnerowicz 
method lOj. By choosing the conformal factor appropriately we can make snre that the 
linearized constraints are satisfied everywhere while also controlling M and J as needed to 
allow ns to prove onr stationarity theorem. 

Onr proof can be snmmarized as follows: Given the spacetime of an axisymmetric perfect 
flnid star that admits a maximal hypersnrface S, we constrnct pertnrbed initial data on S in 
two parts: The first pertnrbation—which we call “the momentnm scaling pertnrbation”— 
scales down the axial part of hij and the polar parts of Kij and Pi while keeping the particle 
nnmber density and entropy density fixed. This is then “fixed np” with a linearized conformal 
transformation so that the linearized constraints are satisfied and so that J, iV, and S remain 
fixed. It is shown that M will be decreased by this pertnrbation nnless the axial part of 
hij and the polar parts of Kij and Pj vanish. The second pertnrbation—which we call “the 
heating pertnrbation”—increases the entropy density in an arbitrarily chosen region in the 
interior of the star while keeping all other variables fixed, and again, the constraints are 
fixed np via a linearized conformal transformation. It is shown that S and M will both be 
increased by this pertnrbation, while N and J remain fixed. Tims, a linear combination of 
these two pertnrbations can be fonnd that increases S at fixed M, N, and J, nnless the 
axial part of hij and the polar parts of Kij and Pj vanish. Tims, if the star is an extremnm 
of S at fixed M, N, and J, these parts mnst vanish. 

Stationarity is then proven as follows. We show that the vanishing of the axial part of hij 
implies the existence of a 0-reflection isometry of S that reverses the direction of the axial 
Killing vector field. Then, the fact that the polar parts of Kij and P* vanish implies that S 
is a “t-0”-reflection symmetric snrface. Bnt it then follows that a neighborhood of S can be 
foliated by (t-0)-reflection symmetric hypersnrfaces, which implies that the neighborhood, 
and therefore the entire spacetime, is stationary. Circularity of the fluid flow follows from 
the vanishing of the polar part of P*. 

In section m we briefly review the Einstein-perfect fluid system. In section m we spell 
out our assumptions on the perfect fluid star spacetimes being considered. In section HVl 
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we show how to decompose the initial data into axial and polar parts, and we write the 
constraint equations as equations on the “manifold of orbits” of the axial Killing vector 
held. In section IVj we construct the perturbation that will be used in our proof, in two 
parts. In section IVTI we state and prove our theorem. Finally, in section I VI It we discuss 
some possible generalizations and extensions. 

Our index notational conventions are as follows: Lower case Latin indices from the early 
alphabet (a, h,c,...) denote abstract spacetime indices, while lower case Latin indices from 
mid-alphabet {i,j, k,...) denote abstract spatial indices. We do not distinguish notationally 
between tensors on a (3-dimensional) spatial hypersurface and tensors on the (2-dimensional) 
manifold of orbits of the axial Killing vector held. Lower case Greek indices ...) denote 
coordinate labels and coordinate components of tensor helds. 

Additionally, three distinct metrics appear in the paper. The (4-dimensional, Lorentzian) 
spacetime metric is denoted gab, and its associated derivative operator is denoted Va. The 
(3-dimensional, Riemannian) spatial metric is denoted hij, its associated derivative operator 
is denoted Di, and its scalar curvature is denoted R{h). The metric on the (2-dimensional) 
manifold of orbits is denoted its associated derivative operator is denoted V^, and its 
scalar curvature is denoted 77(/i). 


II. THE EINSTEIN-PERFECT FLUID FORMALISM 


Here we review the Einstein-perfect huid formalism; for more comprehensive overviews 
see, e.g., { 2 , 11|. The local held variables on a spacetime manifold Ai that describe the 
Einstein-perfect huid system consist of the spacetime metric gab, the particle number density 
n, the entropy per particle s, and the huid 4-velocity satisfying = —1. The stress 
energy tensor has the form 

Tab = ip +p) UaUb + pgab (3) 


where the energy density in the rest frame of the huid, p, is determined by a prescribed 
equation of state. 


P = pin, s). 


( 4 ) 







and the pressure, p, is given by the Gibbs-Duhem relation 


p = —p + pn + Tsn, 


(5) 


where 


is the temperature and 


m 1 

n os 


dp ^ 


( 6 ) 

( 7 ) 


is the chemical potential. The Einstein-perfect fluid held equations consist of Einstein’s 
equation, 

Gab = ^'^Tab, ( 8 ) 


and conservation of the particle number current. 


Va (nu") = 0. (9) 

The latter together with conservation of stress energy (which is, of course, implied by Ein¬ 
stein’s equation) implies conservation of the entropy current, 

Va {snu°‘) = 0. (10) 


The equation of state is required to be such that for all allowed (n, s) we have 


P > 0, P > 0, T > 0, 0 < < 1, 


( 11 ) 


where 




dp 

dp 


dp/dn 


( 12 ) 


dp/dn 

is the square of the speed of sound in the huid. Under these conditions, the Einstein-perfect 


huid system is well posed; see e.g.. 


12| . Additional assumptions on the second derivatives of 


the equation of state should hold so that the huid is locally thermodynamically stable—see 
—but they don’t play a role here. 


As we are interested in conhgurations that represent stars, we will be considering only 
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asymptotically flat solutions in which n has compact spatial support. The fundamental 
conserved global quantities are then taken to be the ADM mass, M, the ADM angular 
momentum, J, the total number of particles. 


N = 



nu°‘h'a^/hd^x, 


( 13 ) 


and the total entropy, 

S = — J snu^'i'aVhd^x, (14) 

where S is any Cauchy hypersurface, is the unit future-directed normal to S, and \/h d^x is 
the induced volume element on S. By ([9]) and flTOll . N and S are conserved, i.e., independent 
of the choice of S. 


Initial data on a Cauchy hypersurface is given (see, e.g., [l2|) by T = (h^, iCjj, n, s, u*), 
where h^j is the induced metric, Kij is the extrinsic curvature, and u* is the projection of 
tangent to the surface. The quantity u'^i'a appearing in ([9]) and flTOD can be written in terms 
of initial data by noting that gabU°‘u^ = — 1 gives 




= —Vi + 




(15) 


where we have dehned 


= hiiU^u^ 


ij • 


( 16 ) 


The initial value constraint equations are 


D^Kij - DiK^. + STrPi = 0 (17) 

and 

- R{h) + KijK^^ - {K\Y + IGwe = 0, (18) 

where the fluid momentum density. Pc = —h'^'hJ^Tab, is given by 

+ P) \/l + Uu 


( 19 ) 
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and the energy density, e = z/“z/^Ta6, is given by 

e = {p + p) [l + V?) - p. (20) 

Clearly, on a maximal hypersnrface—i.e., one with = 0—the constraints rednce to (fTj) 
and ([2|). 


III. ASSUMPTIONS 


Here, a perfect fluid star spacetime, {Ai, gab,n, s,u°‘), is a globally hyperbolic, topologi¬ 
cally R^, asymptotically flat solntion of the Einstein-perfect flnid held eqnations snch that 
n has compact spatial snpport, and snch that the “physical reasonableness” conditions on 
the eqnation of state, flTT|) . hold. 

We will consider only perfect flnid star spacetimes that obey the following additional 
assnmptions: (i) The spacetime is axisymmetric, meaning that the symmetry gronp contains 
a f/(l) snbgronp with spacelike orbits; eqnivalently, that the spacetime admits a spacelike 
Killing vector held with closed orbits, snch that we also have £^n = £^s = £ij,u°‘ = 
0. (ii) The U{1) symmetry gronp acts trivially, in the sense we will define at the end of 
this section, (iii) The spacetime admits a maximal (iCt = 0) asymptotically flat Canchy 
hypersnrface, E. 

While axisymmetry is a gennine restriction (we will discnss the possibility of relaxing 
it in section lylT)) . assnmption (iii) on the other hand is thonght to exclnde only “nnphys- 


ical” spacetimes containing white hole regions, as in Brill’s explicit example. 


13|, of an 


asymptotically flat spacetime that does not admit a maximal hypersnrface. (His example 
consists a “pocket” of expanding dnst glned into the left wedge of the maximally extended 
Schwarzschild spacetime.) In fact, Bartnik [l^ has proven the existence of maximal hy- 
persnrfaces in asymptotically flat spacetimes that satisfy a “nniformity condition” in the 
interior ■ 

For asymptotically flat solntions of Einstein’s eqnation with matter that satisfies the 
strong energy condition, the following two related results are known: (a) If the spacetime 
admits a maximal hypersnrface E, then there exists a unique foliation by maximal hypersur- 


^ Existence of maximal hypersurfaces has been proven for asymptotically flat stationary spacetimes—i 


see 


15 


16l| —but obviously we do not want to assume stationarity here. 
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faces in a neighborhood of S; see [^. (b) If the spacetime is axisymmetric, the axial Killing 
vector held, is tangent to any maximal hypersurface S. This follows from axisymme- 
try, asymptotic hatness, and the aforementioned uniqueness of the maximal foliation: If 0“ 
was not tangent to each maximal hypersurface, one would obtain an inequivalent maximal 
foliation by performing a rotatiorjf . Since the perfect huid matter under consideration here 
satishes the strong energy condition—i.e., {p + p) >0 and {p + 3p) > 0, which follow from 
dnD — assumptions (i) and (iii) imply that we have a unique maximal foliation in a neigh¬ 
borhood of S, and that the axial Killing vector held 0“ is tangent to each hypersurface in 
this foliation. 

As 0“ is tangent to S, we will henceforth denote it as a spatial vector, 0®, on E. Let 
A G T, denote the rotation axis, i.e., the set of points in S at which 0® vanishes, and let 
E = S \ ^. Our assumption, (ii), of trivial action of the U{1) symmetry group is that we 
can write 

E = C> X t/(l), (21) 

where (9 is a 2-dimensional manifold, called the manifold of orbits of Written this way, E 
has the structure of a trivial principal hber bundle with hber group U{1) and base manifold 
O. We believe that under our assumptions that space is topologically and asymptotically 
hat, one automatically has E = x f/(l). We have not attempted to prove this, but if true 
it would make assumption (ii) redundant. 


IV. AXIAL/POLAR DECOMPOSITION, THE MANIFOLD OF ORBITS, AND 
THE CONSTRAINTS 


Many of the constructions in this section may be viewed as a specialization of the “mani¬ 
fold of orbits” formalism of Geroch to the case of a Riemannian manifold. Our exposition 
and notation closely follows js], where the case of axisymmetric vacuum spacetimes in higher 
dimensions is considered. 


On the maximal hypersurface E, we can uniquely decompose the huid momentum density 


® A precise version of this argument was given in the vacuum black hole case in [Sj , which is easily generalized 
to the case where matter satisfying the strong energy condition is present. 
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into its parts parallel and perpendicular to the axial Killing vector field, 0*, via 


Pi — 'P4>i Piy 


(22) 


where (ffpi = 0. Similarly, the extrinsic curvature can be decomposed via 


^ij “t“ ^ij") 


(23) 


where kij and /C* are orthogonal to 0*. We can decompose any axisymmetric tensor field 


nl]- 

kl 


on S into it’s “axial” and “polar” parts via such a decomposition, where the axial 


part consists of any terms in the decomposition that contain an odd number of copies of 0*, 
and the polar part consists of any terms in the decomposition that contain an even number 
of copies of 0*. Thus, the axial part of the fluid momentum density is while the polar 
part is Pi- Similarly, the axial part of the extrinsic curvature is 2/C(j0j), while the polar part 
is kij + Under the action of a 0-refiection isometry—if one exists—the axial part of 

an axisymmetric tensor will reverse sign, while the polar part will remain unchanged. 

As explained in Sj, any axisymmetric tensor field on S that is orthogonal in each of its 
indices to 0* can be projected to the manifold or orbits, O. The tensor fields kij, k, fCi, pi, 
and V are examples of tensors that project to O. As mentioned in section [H we will not 
distinguish notationally between such tensors on S and their projection to O. We would like 
to write the initial value constraint equations, ([T]) and ([2]), completely in terms of quantities 
that project to O. To do this, we need, in addition to the decompositions (12^ and (123|1 
of Pi and Kij, a way to write the spatial metric, and the quantities in the constraints 
associated with it—i.e., Di and R{h )—in terms of objects defined on O. Following Sj, we 
can encode all of the information contained in hij in a triplet, (<F, Aj,/i^), of tensor fields 
that project to O. We now define these quantities. 

First, $ is defined as the norm-squared of 0*, 


$ = 0 , 0 *, 


(24) 


which, as it vanishes in S only on A, is positive everywhere on S. Next, is defined as 


Pij — hij 


(25) 
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which, when projected to (9, is a (Riemannian) metric on O. Finally, to define Rj, we first 
consider the object As explained in Q, $“^0* defines a connection on the principal 

fiber bundle S = (9 x f/(l). Choose a smooth global cross-section, iS, of S, and let A, be the 
axisymmetric 1-form field on E whose pullback to S is equal to the pullback of $“^0* to S 
and is such that = 0. Then Ai contains all the nontrivial information contained in the 
connection but is “gauge-dependent” in that it depends on a choice of cross-section 

S. Given local coordinates on S, one can locally construct coordinates {ip,x^,x‘^) 

on S by using {x^,x^) to label the orbits of 0* and by choosing {d/dtpy = 0*. In these 
coordinates, the spatial metric, hij, takes the form 

= ^dip'^ + 2^Afj,d^pdx^ + + <FA^Ay) dx^dx’^. (26) 


Thus, to show that hij admits a 0-refiection isometry, we must show that the gauge freedom 
in choosing S (equivalently, a coordinate redefinition of the form ip ^ ip + f{x'^)) can be 
used to set 0, so that hij has a “block diagonal” form. 

The curvature 2-form, Fij, is defined as the curvature of the connection i.e., 

Fij = 2D[i (fjjj . (27) 


It is easily verified that Fij is orthogonal to 0*, and therefore projects to O. As shown in 
8|, it follows that 

Fij = 2DiiAj^. (28) 

The Frobenius condition for local 2-surface orthogonality of 0* is that the “twist” of 0* 
vanishes: 

(j)[iDj(j)k] = 0. (29) 

But from fl27l) we have 


4'[iDj4'h] = (44 Vi,|) = 


( 30 ) 


SO the Frobenius condition is equivalent to Fij = 0. 


As shown in 


a 


, both Aj and Fij can be smoothly extended to the rotation axis. A, and 
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therefore the relation 


Fij — 2D[iAj-^ 


( 31 ) 


holds everywhere on S. By assumption, S is simply connected, and hence if F^j vanishes 
then there exists a smooth function y such that Ai = DiX- In that case we can set A* —)■ 0 
by displacing S by —y along the orbits of 0*; equivalently, we can redehne cp ^ cp + x and 
thereby set = 0 in fl26l) . Thus, F^j = 0 is necessary and sufficient for the existence of a 
0-reflection isometry of hij. 

We now have the ingredients in place to write the constraint equations, (1) and ([2]), as 
equations that project to the manifold of orbits. The derivation differs from 8| only by the 
addition of contributions from the fluid matter and in the specialization to a single 0*, so 
we will not repeat the derivation here. The axial part (i.e., the component along 0*) of the 
momentum constraint is 


\/$ 


+ 87r$P = 0, 


(32) 


where is the derivative operator associated with /ijj. The polar part of the momentum 
constraint is 

- ^JC^Fij + 8npi = 0, (33) 

and the Hamiltonian constraint is 


(Vi^) (p*<h)+2<h/cw,+ 


24)2 








+167r£ = 0, (34) 


where 7^(p) is the scalar curvature of pij. 


V. THE PERTURBATION 

Perturbed initial data on E that generates a solution to the linearized Einstein-perfect 
fluid equations on A4 consists of the collection of tensor helds 

54/ = (^6hij,6Kij,6n,6s,6u’‘) (35) 

on E, subject only to the condition that these quantities satisfy the linearized constraint 
equations, i.e., the linearization of flT7)l and flT8|l . 
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We will write our perturbation of interest, as a linear combination of two perturba¬ 
tions, (“the momentum scaling perturbation”) and (“the heating perturbation”), 

each of which satishes the linearized constraint equations individually. These two perturba¬ 
tions can be given the following interpretation: scales down the “non-circular” (i.e., 

polar) part of the fluid motion, thereby decreasing the total energy, M. uses this energy 
to “heat up” an arbitrarily chosen portion of the star, thereby increasing S. Neither of these 
perturbations change J or N. 

Each of these two perturbations will in turn be constructed in two parts: a hrst part which 
satisfies the momentum constraint but not the Hamiltonian constraint, and a second part 
consisting on a linearized conformal transformation designed to “£x up” the Hamiltonian 
constraint. 

The first part of will be designed to scale down the quantities {Fij, , n, pi) by a 
constant factor while keeping the quantities s) fixed, where 

= Vhn\/l + V?. (36) 

(Note: It is not immediately obvious that perturbed initial data, fl3^ . can be chosen to 
accomplish this; we will show that it can.) By inspection, this keeps unchanged both the axial 
and polar parts of the momentum constraint, (|3^ and (l33il . and, since fh and sTl are the 
integrands of flT^ and flTTl) . it keeps N and S fixed. As we will show, it also does not change 
M or J. Furthermore, it does not change the first 4 terms of the Hamiltonian constraint, (l34)l . 
but it decreases the term in square brackets and, as we will show, it decreases the last term, 
IGvre. Thus, the first part of will not satisfy the linearized Hamiltonian constraint. The 
second part of 5 ms' h will fix this—it will be a linearized conformal transformation chosen to 
make the linearized Hamiltonian constraint hold while keeping J, N, and S fixed. We will 
show that this necessarily decreases M unless each of the quantities {Fij, , k, pi) vanish 
everywhere on E. 

The first part of will be chosen to increase s in an arbitrarily chosen region in the 
interior of the star (while keeping it fixed everywhere else) while keeping the quantities 
{hij, Kij, Pi, Dl) hxed everywhere (again, we will show that this can be done). This satisfies 
the linearized momentum constraint but not the linearized Hamiltonian constraint, since, 
as we will show, it increases e in the region where the entropy density is increased. Again, 
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we will fix this up with a linearized conformal perturbation, and the resulting perturbation 
will be shown to increase M and S while keeping J and N hxed. 

Thus, unless the quantities (Fj^, k, pi) vanish everywhere on S, we will be able to hnd 
a linear combination of and that increases S while keeping M, J, and N fixed. 

A. The “Momentum Scaling” Perturbation 

We specify (5 ms^ as a sum of two parts: 


(-'ms ^ (-'ms ^ ' "ms ^ • 


(37) 


(5Ms5k will be the perturbation that scales down {Fij, , k, pi), while (5Ms5k will be the lin¬ 
earized conformal perturbation. We will now specify these perturbations, starting with 

First we specify the perturbed gravitational variables, (5MsFjj), via 


= -kij - K,4>i4>j - 2<FK0(iAj) - 2$/C(iAj), 


or equivalently, in terms of the variables introduced in section IIVI 


= 0, (5«/x,,- = 0, (5« A, = -A, 


= n 


(39) 


Note that we are considering 0* to be hxed as a vector held, i.e., it is unchanged under any 
perturbation we consider here. In particular, (5Ms<?5i = (5 ms(^p(?5^) = = —*hAj, which 

accounts for the last two terms in (j38il . It follows from (j3Til that 


A(1) p.. = —p.. 
"ms-^ o- 


(40) 


Next, we would like to specify the perturbed huid variables, (^msU, (5Mk, ^MsM*) , SO that 
the following three conditions hold: (i) the perturbed huid momentum density satishes 


&i = -p^ - 


(41) 
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or equivalently, in terms of the variables introduced in the previous section, 


■SSb = -B. = 0; 


(42) 


(ii) the perturbed particle number density on S is zero, i.e.. 


= i5S (\/fcnvTTt?) = (nvTTt?) = 0. 


( 43 ) 


where the last equality holds because (|38|1 gives ^msV^ = 0; and (iii) the perturbed entropy 
density on S is zero, i.e., 

= (44) 

which, given (ii), is equivalent to (5 msS = 0. These three conditions uniquely determine 
<5 msS, (5msW*) , and formulas for these can be calculated explicitly from the dehnition of 
Pi, (IT^ - and from the formula for p, ([5]), which in particular implies that 


d f p + p\ I dp 2PPP 

Co 


dn \ n ) ndn ^ 


( 45 ) 


where was defined in (IT^ . We hnd that these formulas 


are 


Cs = 0 


nfpi 


n^^^UiUj 


= - - 

\l + {I - cl)uP^{l + u^){p + PY [l + (1 - c2)M2j 


(46) 




2 i , (-/ + AjYY) 

Ciu'- + Jr -r ; 


n 


(p + p)-\/l + 




n 


ci5P - p% + 71,0* 


uY 


where we have used p* = (p + p)\/l + u'^pPu^ to write these in a way that explicitly shows 
that these dehne a smooth perturbation (i.e., that there is no problem near the edge of the 
star where (p + p) —)■ 0). From (|2n|l we can calculate the perturbation to £ under and 

we hnd 

= -ip + ( 47 ) 


which is manifestly non-positive. 
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By inspection, (5 ms^ prodnces no change in the momentnm constraints, fl3^ and (I33l), 


and prodnces no change in the hrst fonr terms of the Hamiltonian constraint, (jMll . However 


it produces changes in the other terms, giving 


C i-Rih) + + Idvre) = -Q, 


(48) 


where 

p = •F*-’ + 2k^^k-■ 4- H__ ('49') 

Clearly, Q is nonnegative, and is zero if and only if = kij = k = pi = 0. 


Next we specify the linearized conformal perturbation, Our construction follows 

8j, but differs in that we must also specify the perturbed fluid variables and in that we do 


not have a black hole horizon present. The perturbed gravitational variables are given by 




( 50 ) 


where the scalar function will be specihed below. The perturbed fluid variables are 


specihed by requiring 



5<S91 = 0 


( 51 ) 




0 . 


The (unique) values of 5 msS, that give (EH) are 


"MS'’ 


0 


_^ + _ 

2[l + (l-c>2] 


ipn 




2[l+(l-c>=] 


( 52 ) 


and from these, one hnds 

4“^ = -^(3 + 4'U^)(P + P)</’- 


(53) 
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(5ms^ satisfies the momentum constraint, 

^ ,5^=^ [VhD^K.j + STrVhP?) = 0, (54) 

but not the Hamiltonian constraint: 


^MS {-Rih) + + 167r£) 


+ R{h)^p - + 167r(5(,"^e 

2D‘^%jj — 2KijK''^'il) — Svr [(p + 3p) + 2u^{p + p)] 


(55) 


We choose ip to be the solution of 


- D'^tp + [KijK"'^ + 47r(p + 3p) + 9,'Ku^{p + p)] p) 



(56) 


with the boundary condition -0 —)■ 0 at inhnity. As t 


tiplying p) is nonnegative, by standard arguments 18 


le quantity in square brackets mul- 


19| there exists a unique solution 


of (]56|l satisfying this boundary condition. It immediately follows that the perturbation 
(5 ms4/ = <5ms'I^ + (5ms satishes the linearized constraint equations. 


Lemma 1. The perturbation (5ms' h (see fl38ll . fl46l) . (1501) . ((52]), fl56l) ) has 5 msN = 5msS = 
^MsJ = 0. Furthermore, SmsM < 0 with equality holding if and only if Q = 0 everywhere on 

E. 


Proof: The total particle number, f[T3|l . and the total entropy, flTT|) . are 


iv = y 



(57) 


From (1431) and (144j) we have SlffN = ^ms^" = 0. From (|5T1) we have SmIN = 5mIS = 0. Thus, 

6msN = (5ms5 = 0. 

The ADM angular momentum is 


J = f KijpVy/^d'^x, (58) 

8 ^ Js^ 

where Soo denotes a sphere with radius taken to inhnity, r® is its outward pointing unit 
normal, and dfx is its induced volume element. Only the axial part of Kij contributes 
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to this integral, and this is unchanged under (5 ms^ (see W), so we have 5mIJ 
-0 —)• 0 at infinity we also have ^ms'^ = 0. Thus, = 0. 

The perturbed ADM mass is given by 

5M = ^ {Dkdhij - DMjk) ^/7 d^x. (59) 

^ Soo 

Since is purely axial, we have S^lM = oQ, so, substituting we find 

6^sM = S^2M = -^ [ (60) 

Soo 

To evaluate this, let ^ be the solution of 


10 , 


= 0 n. Since 


- + [KijK^^ + 47r(p + 3p) + 87iu^{p + p)] ^ = 0 (61) 


with boundary condition ^ ^ 1 at infinity. By the strong maximum principle 20| ^ is 


strictly positive. Furthermore, by (l56ll . we have 


= (62) 

Integrating this equation over S and using that "0 —)■ 0 and 1 at infinity, we find 

f {r"Di'ilj)y/^(fx = ^ j iQ\/hd^x, (63) 

J Soo ^ S 

so that we have 

f iQ\/hd^x. (64) 

16vr Js 

Since ^ > 0 and Q > 0, we have < 0 with equality holding if and only if Q = 0 

identically. □ 


® Another way to see this is to note that under the assumption of axisymmetry we may write 
J = /j, Vhfjj^Pi = '/h^V, thus, since J^g-v/zi = 5 ms$ = ^msF = 0, we have ^ms J = 0. 

^ One might question whether it makes sense to talk about the ADM mass or angular momentum of a 
perturbation that doesn’t satisfy the constraints. For convenience, we take the attitude here that 6M and 
6J are defined by the formulas (l59l) and (l58l) . independent of whether the constraints are satisfied. For 
our final perturbations 5 ms' 1' and 5 h'I' that do satisfy the constraints, the values given by these formulas 
are then equal to the ADM values. 
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B. The “Heating” Perturbation 


As we did with 5 ms d', we specify as a sum of two parts: 




(65) 


will increase s in an arbitrarily chosen region in the interior of the star while keeping 
Tl and Pi unchanged, and it will not change {hij,Kij). Again, will be a linearized 

conformal perturbation chosen to make the perturbation satisfy the linearized Hamiltonian 
constraint while not changing J. We will now specify these perturbations, starting with 


Firstly, we specify 


S^^^hij = 0 

S'n'Kii = 0 . 


( 66 ) 


Next, we specify the perturbed fluid variables 5 h^s, 5h'w*) by the conditions 


= / 
5 ^'^ 01 = 0 
= 0 , 


(67) 


where / is an arbitrarily chosen smooth function on S that is positive in some arbitrarily 
chosen region in the interior of the star and zero everywhere else. The unique values of 
(5 h ^n, 5 h that reproduce fl^ are 




S"'u‘ 


nv? (riT + g) / 

(p + p)[l + (1 -C2)'il2] 

(l + «^) («T+|)/ 

(p + p)[i + (1 - c'i)u^] 


( 68 ) 


■u 


Note that since we have chosen / to be supported only in the interior of the star where 
(p + p) > 0, these define a smooth perturbation. We can calculate from these formulas, 
and we find 


= riTf, 


(69) 
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which is positive, by our assumption of positive temperature (ITT]) . 

Clearly, satisfies the linearized momentum constraint, but it violates the linearized 
Hamiltonian constraint, giving 


{-R{h) + KijK^^ + IQ-ke) = = 167rnT/. (70) 


Now we specify the linearized conformal perturbation, The perturbed gravitational 

variables are given by 


= i/jhij 


(71) 


where ^|J is chosen to be the solution of 


— D‘^%jj + [KijK^^ + 47r(p + 3p) + 87ru‘^{p + p)]'ijj = SiinT f (72) 

with boundary condition -0 —)■ 0 at infinity. The perturbed fluid variables, (5 h'u, (5h^s, 
are chosen to make (|5T|l hold with ^ which again gives (|5^ and (|53|1 with ^ 
and -0 —t '0- 

With these choices, it follows that (fnlh = satisfies the linearized constraint 

equations. 

Lemma 2. The perturbation (5h4/ has = 5^J = 0. Furthermore, it has 6^3 > 0 and 
5^M > 0 . 

Proof: We have (5 h^T1 = (5h^T 1 = 0, so ^niV = 0. We have = 0 and -0 —t 0 at infinity, 

so from fl58ll we have 5^J = 0. We have (5 h^s = 0, so 

= = J fmd^x>0. (73) 

Finally, we have S^Rij = 0, so = 0, and by the same argument given in the proof of 

lemma [H but with the replacement (—Q) —)• IbvrnT/, we have 

6^M = = [ ^nTf^/hd^x > 0. (74) 

Jt. 

□ 
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VI. STATIONARITY THEOREM 


Theorem. Let {Ai, gab,n, s,u°‘) be a perfect fluid star spacetime satisfying the assumptions 
of section \II^ with axial Killing vector field Suppose that 5S = 0 for all solutions 

of the linearized Einstein-perfect fluid eguations that have 6M = 6J = 6N = 0. Then 
{At, Pabyn, s,u°‘) is stationary with circular flow — i.e., there exists an asymptotically timelike 
Killing vector field such that £t(f)°‘ = £tn = £ts = £tu‘^ = 0, and such that is a linear 
combination of K and 

Proof: Let S be a maximal Cauchy surface, and consider a perturbation having the initial 
data 

5^ = 5h^ + A5ms^, (75) 

where cind were constructed in the previous section, and A is a constant. If 

Q is nonzero anywhere on S then by lemma [1] we have 6msM < 0, so we can choose 
A = 5 h^/|^ms^| which makes dM = 0. But by lemmas [T] and [2] we have dN = dJ = 0 and 
> 0, which is a contradiction. Thus, Q = 0 everywhere on E, where Q was dehned by 
(09]). Consequently, we have 

Fij = kij = K = pi = 0 (76) 


everywhere on E. 

As shown in section IIVI Fij = 0 implies that there exists a (/)-reflection isometry is : 
E —)■ E that reverses the sign of the axial Killing vector held (i.e., i^{4>fl = —0*)- Then 
kij = K = Pi = 0 implies that Kij and Pi are odd under is (be., i'^{Kij) = —Kij and 
i^{Pi) = —Pi)- It follows from the dehnition of Pi, flT^ . that K is odd under is. 

Now consider the “f-0”-rehection, i : M —>■ M, obtained by mapping a point p lying at a 
proper time r along a normal geodesic starting at s G E to the point q lying at proper time 
—r along a normal geodesic starting at is(s). Since i maps the initial data (hij, Kij,n, s, K) 
on E into itself, by uniqueness of Cauchy evolution, i is a symmetry of the solution. 

Next consider the foliation by maximal hypersurfaces E(t) in a neighborhood of E = E(0); 
by such a foliation is guaranteed to exist and is unique. By repeating the above argument 
for each surface in the foliation, we have a (t-0)-rehection symmetry, fl, about each E(f). 
We can choose the parameter t so thaio h[E(t + s)] = E(t — s). The desired 1-parameter 


The uniqueness of the maximal foliation allows us to make this choice. To see this explicitly, one can 
choose t to be the proper time at which a particular normal geodesic starting at S passes through S(t). 
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group of time-translation symmetries, 4 ; —?• A4, is then obtained via composition of 

reflections: 

It = it/2 o ^0 (77) 


(it can easily be seen that It takes S(s) to T,{t + s)). We have thus shown stationarity in a 
neighborhood of S, but, by uniqueness of Cauchy evolution, this implies stationarity on all 
of Af. 

Finally, circularity of the fluid flow is shown by noting that the fact that u* is odd under 
is means that u* is parallel to (/>*. Thus, is a linear combination of 0“ and the unit normal 
vector held, 12 °', to S(t). But must also be a linear combination of z/“ and 0“ (having a 
component tangent to S and orthogonal to 0* would contradict the fact that is odd under 
i), thus is a linear combination of and 0“. □ 

The following corollary was proven in j^j®: 


Corollary. Under the hypothesis of the above theorem, the quantities hi, T = |17|T, and 
ft = \V\ fi are uniform, where 




(78) 


In other words, axisymmetric thermal equilibrium stars are rigidly rotating with uniform 
redshifted temperature and uniform redshifted chemical potential. Furthermore, the hrst law 
of thermodynamics, 

6M = fdS + ftdN + nSJ, (79) 


holds for all solutions of the linearized Einstein-perfect fluid equations. 


VII. GENERALIZATIONS AND EXTENSIONS 

An obvious generalization we would like to be able to make is to relax the assumption 
of axisymmetry in our stationarity theorem. However, having an axial Killing vector held 
was essential in our method of proof, since without it we don’t know how to decompose 
the huid how into a part that does not contribute to the angular momentum (polar) and a 
part that does contribute (axial), so that we can construct a perturbation that slows down 

® The uniformity of ti, T, and fi was shown earlier by 
22 |. 


21j . and the first law then follows from a result of 
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the fluid flow without changing the angular momentum. One might therefore try to argue 
directly that there are no non-axisymmetric extrema of entropy. At first glance, the following 
argument would seem to do this, but, as we will explain, it does not actually work: 

For an axisymmetric star, perturbations consisting purely of gravitational radiation can¬ 
not change J to first order (see, e.g., the discussion surrounding equation (63) in ). But for 
a non-axisymmetric star, they can. Thus, one might think that for a non-axisymmetric rotat¬ 
ing star, one could always construct a perturbation that has 55' > 0 with 5M = 6J = 6N = 0 
by: (i) slowing the rotation, (ii) using most of the energy thereby gained to heat up the star, 
and (iii) putting the rest of the gained energy and all of the gained angular momentum 
into a gravitational wave-packet very far away from the star. By putting the wave-packet 
arbitrarily far away, it can carry the given amount of angular momentum while having arbi¬ 
trarily small energy, thereby guaranteeing that there will be enough energy left over to heat 
up the star. 


However, it has been shown 23| that there exist asymptotically flat non-axisymmetric 


solutions of the initial value constraints that are exactly axisymmetric outside of a compact 
region. For such solutions, a wave packet that is too far away from the star would be in 
the axisymmetric region and therefore could not change J to first order. Thus, the above 
argument does not work, and currently we do not see how to proceed without the assumption 
of axisymmetry. 

Other generalizations include considering (i) spacetimes of dimension d > 4 that are 
axisymmetric—meaning that if there are N < nonzero independent angular mo¬ 

menta (J^,..., J'^), then there exist N mutually commuting axial Killing vector fields 
..., associated with them—and (ii) asymptotically anti-de Sitter spacetimes. As 
shown in [sj, as long as the action of the rotational symmetry group is trivial in the sense 
discussed at the end of section im (in the higher dimensional case this is believed to be a 
genuine restriction, as discussed there), one can construct a smooth perturbation on S that 
scales the curvature 2-form and the polar parts of the extrinsic curvature, where now these 
are objects that have extra “internal” indices corresponding the the multi-dimensional space 
of rotational symmetries. Thus, the constructions of section |V] should generalize straight¬ 
forwardly to the case of higher dimensional stars having trivial action of the rotational 
isometries. Furthermore, as noted in js], including a negative cosmological constant does 
not affect our ability to solve the linearized Hamiltonian constraint fl56|) of a conformal per- 
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turbation with a given source, so our theorem should generalize straightforwardly to stars 
in asymptotically AdS spacetimes. 

Finally, we discuss the analog of the conjecture of section |T] for the case of black hole 
spacetimes. Namely, that if the area of a cross section of the event horizon—which plays the 
role of entropy in black hole thermodynamics—is an extremum at hxed M and J, then the 
black hole is stationary and axisymmetrico. If the area of a cross section of the event horizon 
is an extremum, then the expansion of the null geodesic generators must vanish at that 
cross section, because otherwise an inhnitesimal diffeomorphism along the null generators 
(i.e., a “time translation”) would increase the area, contradicting it being an extremum. 
But by Hawking’s area theorem 2^, the expansion must be nonnegative everywhere, so 
the Raychaudhuri equation implies that the expansion and shear vanish everywhere on the 
horizon (in particular, this means the area is constant). It might then follow by arguments 
similar to those used to prove the rigidity theorem 2J, [2^ that there exists a Killing vector 
held normal to the horizon, but we don’t have a proof of this. 

The above argument—which uses the “global methods” of the area and rigidity theorems— 
obviously has a very different character than the “scaling of initial data” method we used 
to prove our huid star stationarity theorem. This can be understood by the fact that 
the event horizon is a globally dehned quantity and cannot easily be located from initial 
data. However, apparent horizons are dehned in terms of initial data on a hypersurface, so 
one might try use our methods to prove stationarity of axisymmetric extrema of apparent 
horizon area on a maximal hypersurface with a trapped region. As the apparent horizon 
and event horizon coincide in the stationary case, this seems like a reasonable conjecture to 
make. But to give it the interpretation of “extrema of black hole entropy are stationary,” 


one would have to argue that apparent horizon area, rather than event horizon area, s 
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lould 
, but 


play the role of black hole entropy. Such arguments have been made, e.g., in [26 
we will not take a position here. In any case, it appears that it is not actually possible to 
prove such a conjecture using our methods for the following reason: The momentum scaling 
and conformal perturbations move the location of the apparent horizon within the maximal 
slice, and “moving it back” (via an inhnitesimal diffeomorphism) generally changes the area 
in an uncontrolled way. 


Note that the converse—i.e., stationary black holes are extrema of event horizon area at fixed M and 
J —follows directly from the first law of black hole mechanics. 
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